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1. Introduction 



e + e~ annihilation into hadrons proves to be a wonderful laboratory for detailed experi- 
mental tests of QCD. The high statistics data collected from hadronic Z° decays in LEP-1 
and SLD experiments |]] allow one to perform detailed studies of perturbative QCD and 
to reduce the domain of our ignorance on the physics of confinement. The data have 
convincingly demonstrated the dominant role of the perturbative phase of jet evolution 
and supported the hypothesis of local parton-hadron duality ]2|J3[]. We have now quite 
successfully entered the stage of quantitative tests of the so-called Modified Leading Log 
Approximation (MLLA) which allows one to calculate systematically the inclusive jet 
characteristics up to terms of relative order ^fal In some special cases next-to- 

MLLA corrections are also calculated, (for reviews see ||||[7]]). 

In November 1995 LEP was operated at centre-of-mass energies 130-140 GeV (LEP- 
1.5). The first experimental data [§,|yn|[ll] show good agreement with perturbative 
expectations. In particular, semihard QCD results have become available now. 

Further studies of the various aspects of QCD dynamics will be provided by LEP-2 
and future linear e + e~ colliders (for discussion, see e.g. [l2j and [|L|) as well as by 
HERA and TEVATRON. It is worthwhile to mention the unique opportunity to compare 
the particle spectra measured by the same experimental group at the same accelerator at 
different energies. 

The aim of this paper is twofold. First of all we would like to perform the complete 
analysis of the recent data on inclusive charged particle distributions collected by LEP- 
1.5. We show that the inclusive charged particle spectra and mean multiplicity agree 
impressively well with the analytical QCD results. Secondly we would like to discuss the 
working tools relevant for the detailed treatment of the forthcoming data. 



2. Inclusive charged particle spectrum in QCD jets 

In the MLLA this spectrum can be obtained as solution of the appropriate evolution 
equation ||f4| in terms of two parameters, the QCD scale A and the k± cut-off Q in the 
cascade. In the case when both parameters coincide (Qo = A) one obtains the so-called 
limiting parton spectrum ||6||l4f which has been proven to be very successful in fitting the 
experimental data of charge particle production in QCD jets. 

Until now, characteristics of an individual quark jet have been studied in greatest 
detail experimentally in the reaction e + e~ — ► hadrons. An inclusive momentum spectrum 
here is the sum of two g-jet distributions. In terms of the limiting spectrum one obtains 
for the distribution in the variable £ = log 1/x with x = 2Eh / \/i 

= 2K h D^,Y) (1) 

where K h is the hadronization constant, the total cms energy and Y = log(-y/s/2A). 
The limiting spectrum is readily given using an integral representation for the confluent 
hypergeometric function ||,|[]: 
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Here a = oto + i£ and cto is determined by tanh ao = 2£ — 1 with £ = 1 — 4. 

is the modified Bessel function of order 5, where B = a/b, a = lL/Vc/3 + 2nf/3Nc, 

6 = (HJVo - 2n/)/3, with n/ the number of flavours and C F = (iVg - l)/2iV c = 4/3. 

The analysis of charged particle spectra using this distribution [|i5j|il||T7ll yields values 
for the effective scale parameter A = A c ^ in the range A c ^ ~ 250 -j- 270 MeV. If both 
parameters Qq and A are kept free in the fit one recovers the limiting spectrum with 
Q = A as best solution JET). 

It proves to be very convenient (see e.g. [ ,[T7|] ) to analyze inclusive particle spectra 
in terms of the normalized moments 
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where Me is the mean multiplicity in the jet, the integral of the spectrum. These mo- 
ments characterize the shape of the distribution and are independent of normalization 
uncertainties. The theoretical predictions for the moments from the Limiting Spectrum 
are determined by only one free parameter A c ^. Also one defines the cumulant moments 
K q (Y, A) or the reduced cumulants k q = K q /a q 



19 1 which are related by 
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where the third and fourth reduced cumulant moments are the skewness s and the kurtosis 
k of the distribution. If the higher-order cumulants (q > 2) are sufficiently small, one can 
reconstruct the ^-distribution from the distorted Gaussian formula, see 

ra 



The cumulant moments can be obtained from 
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where 7 w (a s (y)) denotes the anomalous dimension which governs the energy evolution of 
the Laplace transform D U (Y) of the ^-distribution D(£, Y). Equation (§) shows the direct 
dependence of the moments on a s (Y) and thereby implicitly on nj. 

The analytical procedure of |n| allows one to arrive at the following expression for 
the moments: 
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where the parameter B is introduced above and the variable z is given by z = yl6N c Y/b; 

Pq 9 ^ and P^ are polynomials of order 2(q— 1) in z. Explicit results for the full expressions 
for q < 3 can be found in |14| and for q = 3,4 in 



In the calculations, the partons were taken massless and the p t cut-off Q = A ch was 
introduced for regularization; on the other hand, the observable hadrons are massive. One 
can make the simple assumption that the cut-off Qo is related to the masses of hadrons. As 
a first stage in the discussion of charged particle spectra, one can relate Q to an effective 
hadron mass, then for both partons and hadrons E > Q . A consistent kinematical 



behaviour of parton and hadron spectra can be obtained through the relation [[17 



dn(£ E ) = dn{£ E ) 
dp h P dp p 



Eh — — — = E p — — — ; E h = E p > Q c 



at the same energies E^ or £ E . = log(y / s/2P,) (not momentum!), where E h = ^jp\ + Q\ 
and Ep = p p . With this choice both spectra vanish linearly for E — > Qq or £g — > Y. 



3. Comparison with data 

LEP-1.5 opens a unique opportunity to compare the high energy data collected at the same 
accelerator by the same Collaboration^. Later on such a comparison will be continued at 
LEP-2. 

Momentum spectrum 

Figure la shows experimental data on charged particle inclusive momentum distribu- 
tions, dn/d^ p as a function of £ p = log(\/s/2p), obtained by the OPAL Collaboration at 
cms energies y/s = 91.2 GeV (LEP-1) pf and ^/s = 133 GeV (LEP-1.5) 0. Theoreti- 
cal predictions of the Limiting Spectrum ([|) at the same cms energies are also shown; at 
both cms energies, the A c h parameter has been taken equal to 270 MeV, as suggested by 
the moment analysis performed in |T7| , whereas the free overall normalization factor has 



been fixed to the value K h = 1.31^ In view of the forthcoming run at LEP-2, theoretical 
predictions for ^/s = 200 GeV with the same choice of parameters are also shown. In 
Figure lb data on inclusive momentum distributions extracted at LEP-1.5 cms energy by 
different experimental Collaborations, i.e., ALEPH §, DELPHI §| and OPAL QTT 
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compared to the same theoretical predictions as in Figure la. The Limiting Spectrum 
predictions with this choice of parameters reproduce very well the experimental shape 
around the peak at LEP-1 cms energy; the agreement at LEP-1.5 is reasonable, even 
though deviations are visible in the region around the maximum. However, this effect 



s at lower energies such a comparison has been performed by TASSO Collaboration at PETRA Q. 

^OPAL Collaboration ||| found a best value of A ch = 253 MeV at LEP-1 with K h = 1.28 from a 
fit to the shape, whereas A c h — 263± 4 MeV was obtained from a fit to the peak position by using the 
asymptotic formula |ll]]. The value for A c h — 270 MeV has been chosen in ]l7| to reproduce the moments 
over a large energy interval. With a restriction to the LEP-1 data, a smaller value of A c h would fit slightly 
better. 

"DELPHI published the charged particle inclusive energy distribution obtained by assigning to all 
particles the pion mass. The inclusive momentum distribution has been extracted from experimental 
data by using the same assignment to the mass of all particles. 
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could be well attributed to statistical and/or systematic uncertainties, as suggested by 
the large spread of experimental points shown in Figure lb. 

Looking first at the large x p (small £ p ) region, the experimental data can be well 
described at both cms energies by the Limiting Spectrum prediction (P, contrary to the 
simple distorted Gaussian parametrization (not shown). The observed agreement of the 
Limiting Spectrum prediction in a kinematical region where the MLLA approach is not 
expected a priori to be valid is due to the fact that the approximate expression for the 
anomalous dimension derived within MLLA at small x p turns out to mimic reasonably 



well the expression valid at large x p (see |2^JT6|j2^| for more details) 



In the small x p (large £ p ) region, data show a tail which is not well reproduced by 
theoretical predictions taking £ p as an argument in (fj). This discrepancy is actually due 
to kinematical effects |p!6[p7| . As previously discussed, theoretical predictions vanish for 



E < Q , i.e., for > Y. This kinematical effect can be taken into account in a simple 
way by using the relation (|j) between parton and hadron spectra. The corresponding 
theoretical predictions for the spectrum can then be written as: 

~^ = 2K^D^ E ,Y) , 6, = log f (9) 



For E ^> Q these differences vanish and dn/d^ p ~ dn/dt^E, ^e — £ P - It is seen in Fig- 
ures la and lb (dashed lines) that they closely follow the large £ tail of the experimental 
data. The best fit here is obtained with K h = 1.34. It is also worth noticing that a 
common normalization factor, K h , has been used at both cms energies; this result is 
consistent with the expectations of LPHD. 
Moments 

The moments < £ 9 > and the cumulants K q are determined from the spectra Edn/dp 
as a function of the average multiplicity Me is obtained as the integral over of 
the full spectrum Edn/dp. For the unmeasured interval near £^ ~ Y (small momenta) a 
contribution was found by linear extrapolation as in Notice that Me coincides with the 
usual particle multiplicity, M, the integral over £ p , at asymptotic energies. At LEP-1 
(LEP-1.5) the difference is about 10% (8%)0. The cumulants up to order q=4 together 
with the corresponding MLLA predictions at cms energy yfs = 133 GeV are shown in 
Table 1. Notice that DELPHI has also presented data at the Z° peak from radiative events 
collected at LEP-1.5. Cumulants extracted from different experiments are consistent 
within experimental uncertainties; it is important to stress that also radiative events at the 
Z° peak collected at LEP-1.5 are completely consistent with LEP-1 results. Errors quoted 
in the Table are statistical errors only. Note that moments of order q > 1 are independent 
of the overall normalization, which is the main source of systematic uncertainties. On 
the other hand, a systematic error of the order of 0.5 should be taken into account for 



the average multiplicity at LEP-1.5 energy. Also shown in the Table are the results [fL7 
obtained at LEP-1 for comparison. Theoretical predictions are shown for the Limiting 
Spectrum with A c h = 270 MeV; for the average multiplicity Me the first number in 
the Table refers to the prediction (Bl) normalized at LEP-1, the second number to the 



extrapolation of the fit Me = CIo^Mls + c 2 to the energy region 3-91 GeV (T^j. The 



*Good fits have been obtained also for the standard charged particle multiplicity, A/" ]^ i p| i p^|JTl| | . 
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agreement between experimental results and theoretical predictions is very satisfactory; 
the deviation visible in the prediction of the average multiplicity at LEP-1.5 is well inside 
the errors if one correctly takes into account systematics. 
Peak position 

An easily accessable characteristic of the ^-distribution is its maximum £* which 
has been extensively studied by the experimental groups (for a recent review, see e.g. 
The high energy behaviour of this quantity for the Limiting Spectrum is predicted 



as 
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with the constant term given by 
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0.2915 (0.3513) for n f = 3(5). 



Alternatively, one can compute the maximum £* from the Distorted Gaussian approxi- 
mation: 

e=l~\sa (11) 

using the full expression (|7|) for the moments defined in @. One neglects here the 
contribution of cumulants of order greater than 4 to the position of the maximum. The 
third possibility is to extract numerically the actual position of the maximum £* from 
the Limiting Spectrum (|2|); the comparison between the three choices can give us some 
information on the validity of the approximations performed in the first two cases. A 
comparison of these three possible approaches is shown in Figure 2a for rif = 3. Notice 
that since we are plotting vs. Y, there is no dependence of theoretical predictions on 
Kh- 

One can conclude from Figure 2a that the energy dependence of the maximum £*, i.e., 
the slope of the curve, is very similar in the three finite shift, however, is present; 

for example the true maximum of the Limiting Spectrum £* can be expressed in terms of 
the approximation ([H]) as: 
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(The numerical term actually decreases slightly from 0.11 at y/s = 10 GeV to 0.086 at 
130 GeV). 

In Figure 2b we compare experimental data on the maximum £* extracted from 
Gaussian or Distorted Gaussian fit^J] to the central region of the inclusive momentum 
spectrum [0j2|j2|j23j2|,|2|g,p,p,|l|,|,0 with the theoretical predictions obtained 
extracting the maximum from the shape of the Limiting Spectrum. We found a good 
agreement for A c h = 270 MeV, which also describes the energy evolution of the moments. 
The best value of A c ^ which fits the experimental data is slightly different, depending 



^We used the £* values from the original publications or from the fitsEj 
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on which theoretical formula is beeing used. For instance, by using the asymptotic for- 
mula (12), a smaller value of A c ^ would have been better. The difference in the A c ^ value 
is of the order of 20-30 MeV. 

Another possible source of ambiguity is the number of active flavours to be used in 
the theoretical formula. This problem can be studied most easily [0 for the moments 
of the ^-distribution whose evolution with energy follows from (||). In this formula the 
number of flavours enters through the running coupling a s (y,nf). The moments evolve 
at low energy with 3 active flavours and with 4 and 5 flavours after passing continuously 
the respective thresholds. As the MLLA is based on the one-loop expression for a s there 
remains a scale ambiguity. The simplest approach would be to put the thresholds at the 
heavy quark masses, i.e., to set the scale of running a s to However, let us remind that 

a s depends on the transverse momentum k t and kinematics forces k t < \^-- This would 
suggest to move the thresholds to 4 itlq or towards even larger values, since the effective 
value of kt would be in general smaller than the maximum allowed value. Taking the 
presence of heavier flavours into account in this way, one finds small deviations at LEP-1 
energies in the higher moments q > 2 |IT|| . However, for £* using the approximation (|IID 
the effect of the higher flavours is only of the order of 1% and taking rif = 3 is a very 
good approximation. This is demonstrated in Table 2 which compares the predictions for 
£* under different assumptions on the flavour composition and effective thresholds. 

There is also a small difference between the maximum £* of the Limiting Spectrum (|2]) 
and the modified distribution @ by about 1%. This dependence may be considered as 
systematic uncertainty of the predictions. 

Let us also point the attention to the fact that the prediction of the Limiting Spec- 
trum ( |TTD describes the data quite well without the need of any additional constant 
contribution which might be due to the hadronization process. 



4. Conclusions 

The analysis of the recent LEP-1. 5 data demonstrates that the analytical perturbative 
description of inclusive particle distributions in QCD jets is in a quite healthy shape. 
The charged particle spectrum (up to the overall normalization) and the moments with 
q > 1 as well as their energy dependence can be described with only one adjustable scale 
parameter A c /j. For these observables, no additional sizeable effects from hadronization 
are visible. The present data allow one to observe even relatively small effects caused 
by higher order terms. Finally, let us emphasize that in our view it is quite impressive 
that even a limited sample of hadronic events (corresponding to integrated luminosity of 
about 5 pb~ l ) allows one to perform so good tests of the perturbative picture of multiple 
hadroproduction in jets. Having in mind the prospects of QCD studies at LEP-2, this 
looks rather promising. 
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6. Figure Caption 



Figure 1. a. Charged particle inclusive momentum distributions at y/s = 91.2 GeV (LEP- 
1) [[OJ (diamonds) and y/s = 133 GeV (LEP-1.5) (Tl[ (triangles) as measured by 
the OPAL Collaboration in comparison with theoretical predictions of the Limiting 
Spectrum with A ch = 270 MeV (solid line). Dashed lines show the predictions of the 
Limiting Spectrum after correction for kinematical effects. Theoretical predictions 
are also shown for y/s = 200 GeV. 

b. Charged particle inclusive momentum distribution at y/s = 133 GeV (LEP-1.5) 
measured by ALEPH f§ (diamonds), DELPHI @ (squares) and OPAL Collabora- 
tions |nj (triangles), and compared to the same theoretical curves as in a. 

Figure 2. a. Maximum of the inclusive momentum distribution as a function of Y = 
log 2A^~! comparison between different theoretical predictions: maximum numeri- 
cally extracted from the shape of the Limiting Spectrum (solid line), asymptotic 
formula ( |I0D (dashed line), eq. ([□]) (dotted line). 

b. Maximum of the inclusive momentum distribution £* as a function of Y = 
log 2A^"i comparison between experimental data (see text) and theoretical predic- 
tion numerically extracted from the shape of the Limiting Spectrum (solid line); A c h 
= 270 MeV. Crosses are put at cms energies y/s = 200 GeV and 500 GeV. 
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7. Table Caption 



Table 1. The average multiplicity Me, the average value the dispersion a 2 , the skewness 
s and the kurtosis k of charged particles' energy spectra Edn/dp as a function of 
with A c h = 270 MeV extracted from experimental data at cms energies yfs = 91.2 
GeV and 133 GeV. In brackets the theoretical predictions of the Limiting Spectrum; 
the second entry in the average multiplicity column contains the results of the two 
parameter formula Me = C\~IMls + C2; the first one the fit with C2 = 0. 

Table 2. Dependence on the number of active flavours of the theoretical prediction (|TTD for 
the maximum of the inclusive momentum spectrum. Results at three different cms 
energies are shown. 
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